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I. INTRODUCTION 

Quantum mechanics predicts that a particle in a peri- 
odic potential will undergo oscillatory motion in response 
to a constant force. This rather counterintuitive phe- 
nomenon, known as Bloch oscillations (BO), was contro- 
versial for a long time, but has now been experimentally 
demonstrated in very clean semiconductor superlattices 1 
and Bose-Einstein condensates^. 

Can Bloch oscillations also exist in magnetic systems? 
Kyriakidis and LosPl discussed this possibility. They 
considered a system where the particle is a propagating 
domain-wall excitation in an easy-axis one-dimensional 
ferromagnet, and concluded that magnetic BO should 
indeed exist. In particular, the blue crystalline material 
C0CI2 • 2H2O was identified as a promising candidate for 
observing magneti c BO . Searches using neutron scatter- 
ing were performed^, but did not find evidence of BO. 

Another very similar system where one can expect 
BO of magnetic domain-walls is the Ising model in a 
magnetic field having both longitudinal and transversal 
components^. Such a model is believed to be realized in 
CoNb2 06 where indeed an intriguing frequency spectrum 
has recently been observed^. However, in the region of 
momentum space where one expects to find the quan- 
tum mechanical spectral signature of BO — a spectrum 
with equidistant energy levels, the so-called Wannier- 
Zeeman ladder (WZL) — the spectral weight in the ex- 
periment in Ref. is dominated by a strong feature at- 
tributed to additional couplings in the Hamiltonian, a 
"kinetic bound state" , stabilized by next-nearest neigh- 
bor interactional. Thus it appears that additional terms 
in the Hamiltonian prevent BO in CoNb 2 06. 

This might also be the case in C0CI2 • 2H 2 where 
the "kinetic bound state" will be generated by a near- 
est neighbor spin flip exchange coupling that, indeed, is 
present in CoCl 2 • 2H 2 C 4 , but neglected in Ref. [3J It is 
the main goal of this article to investigate its influence 
on the WZL in CoCl 2 • 2H 2 0. 

Our results show that the WZL is present in certain 
regions of momentum space also in the presence of the ex- 



change couplings. However, the neutron scattering spec- 
tral weight of the WZL in CoCl 2 • 2H 2 is less than 1% of 
the total spectral weight at these momenta, thus making 
it difficult to observe the WZL at zero temperature in 
inelastic neutron scattering experiments. 

At finite temperatures, the neutron scattering signa- 
tures of the WZL are more pronounced. We find that 
a relatively high temperature is favorable as the number 
of domain-walls performing BO are exponentially sup- 
pressed with temperature below the largest ferromag- 
netic coupling. Unfortunately, a high temperature leads 
also to collisions of domain-walls that destroy BO. For 
C0CI2 • 2H2O, we find that this collision rate is deter- 
mined by the velocity of the "kinetic bound state" and 
the distance between domain-walls bound states that gets 
smaller as the temperature increases. A way to alleviate 
this is to increase the applied magnetic field that leads 
to a reduced collision rate. However, a large magnetic 
field makes the intrinsic signature of BO weaker as it re- 
duces the BO amplitude. In searching for a compromise, 
we find significant neutron scattering signatures of BO 
at finite frequencies in CoCl 2 • 2H 2 at the temperature 
T = J z /2, magnetic field h z = 0.2J Z , and momentum 
transfer q — ir/2 where the first Bloch mode at lub = 2h z 
carries about 12% of the total spectral weight. 

While we focus on the material parameters for C0CI2 • 
2H2O, our results are analytic and can, with minor ef- 
forts, also be used for aiding searches for magnetic BO 
in other similar materials. 



II. HAMILTONIAN 

We start with the spin-1/2 XYZ ferromagnetic Hamil- 
tonian for a chain with nearest-neighbor coupling in a 
magnetic field 

H = — (j x SfSf +1 + JySf Sf +1 + J z S*Sl +1 + h z S*) , 

i 

(1) 
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which can be written as 

H = H Z +H a + H ± , (2) 

where 

h z = -j z ^2 s i s t+i s fi 

i i 

H a = — J a ^ {SfSf +1 + S i S i+l ) , 

i 
i 

with J a — (J x — J y )/4 and J± = (J x + J y )/^, and 
Sf = Sf ± iSf are the usual raising and lowering op- 
erators. We will assume that J z is by far the largest 
coupling thus the system has an easy-axis and its be- 
havior is Ising-like. The ferromagnetic coupling causes 
neighboring spins to align their ^-components, and with 
a large J z , the excitation energy of a state is mainly de- 
pendent on the number of anti-aligned spin neighbors, or 
domain- walls. Each domain wall costs an energy of J 2 /2 
(our J z is twice that of Ref. [9]). Thus the ground state 
is approximately the ferromagnetic state where all spins 
are aligned along the z-axis. In the absence of a mag- 
netic field and other couplings, the first excited state has 
a single domain-wall, see Fig. [T] This domain-wall can 
be placed between any of the spins, implying a macro- 
scopic degeneracy. The 7J -term lifts this degeneracy re- 
sulting in a band dispersion describing the dynamics of 
a single domain-wall. This mode was first predicted by 
Villain 10 and was subse quently observed in neutron scat- 
tering exp eriment i 1 1 * 1 2 \ 

Kyriakidis and Loss treated a single domain- wall in the 
presence of a finite magnetic field and predicted BCP. 
The BO are caused by H a together with the magnetic 
field that causes the domain wall to oscillate. While the 
single domain- wall approximation is presumably good on 
short time-scales where collisions between domain-walls 
can be ignored, it will break down at longer time scales. 
In a finite magnetic field, this time-scale is likely to be 
very short as a domain- wall and an anti-domain-wall are 
closely bound together. The energy cost of a domain 
of spins anti-aligned with the field is proportional to 
the magnetic field times the number of spins in the do- 
main. Thus the magnetic field induces a linear potential 
between a domain-wall and an anti-domain-wall which 
confines them in a bound state. Therefore the low en- 
ergy excitations will not be isolated single domain-walls, 
but rather bound states of domain- wall/antidomain- wall 
pairs that define the boundaries of a spin cluster of over- 
turned spins, see Fig. [2] The far-infrared light absorp- 
tion experiments on the quasi-one-dimensional material 
C0CI2 • 2H2O in a magnetic field hav e bee n explained in 
terms of such spin cluster excitations^^. In the bound 
state picture, H a and the magnetic field cause the bound 
state to shrink and expand. This is what gives rise to the 
BO and the WZL. 



FIG. 1: A spin state with a single domain- wall. 
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FIG. 2: Two bound states. 

Going beyond the single domain-wall approximation 
also allows the inclusion of the spin-flip exchange Hamil- 
tonian H±. The action of H± on the single domain- wall 
state shown in Fig. [T] produces a high energy state hav- 
ing three domain-walls. In contrast when H± acts on 
the minimal bound state — a single overturned spin, see 
Fig. fright — it can move the whole bound state without 
introducing extra domain-walls, thus H± acts directly 
in the low energy Hilbert space of a single bound state. 
This minimal bound state is the analogue to the "kinetic 
bound state" in CoNb 2 06- 

In zero magnetic field, the predictions for the neutron 
scattering dynamic structure factor were shown to be 
independent^ of whether one considered noninteracting 
domain-walls^ or bound states^. However, as argued 
above, this does not hold in a finite magnetic field, at 
least not on longer time scales. In this article we will 
treat the single bound state exactly in the low energy 
subspace to two domain- walls and assume that thermally 
the system can be well approximated by a noninteracting 
gas of such bound states. 

III. QUANTUM MECHANICS OF A SINGLE 
BOUND STATE 

Let us represent a bound state of a domain-wall and an 
anti-domain-wall as the state 

|i,i) = |...ttU...4.tt...) s (3) 

where the index j — 1, 2, ... N gives the starting posi- 
tion of the down-spin cluster and I = 1, . . . N describes 
its length. N is the total number of spins in the chain, 
which we will take to be a macroscopic number. We will 
extend this representation of states to I — in order to 
also include the ferromagnetic state \j, 0) , which is inde- 
pendent of j. The action of the Hamiltonian on such a 
state can be written as 

H\j,l) =(l-6 L0 ){j z \j,l)+h z l\j,l) 
-J a [\j,l + 2) + \j -2,1 + 2) 
+ l-2) + \j + 2,l- 2)) (1 - 5 L2 )(1 - 8 lA )] } 
-^(1^2)^,0 + ^0)^,2) 

-Jx(\j + 1,1) + |j - 1, (4) 

where we have neglected terms that result in more than 
two domain- walls. 
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FIG. 3: Examples of the actions of H a on particular states 
from the and 2 domain- wall sectors. Only cases that yield 
the minimal amount of increase in domain walls are shown. 
Operating with H a are shown on the left, while the right hand 
side shows the effects of operating with H±. The increase in 
the number of domain walls is indicated to the right of each 
process. Note that H a has the ability to move domain walls 
without increasing their number when acting on a state with 
one or more domain walls (left side, the two lowest processes), 
while H± lacks this ability with the exception that it can move 
a single overturned spin without creating new domain walls 
(right side, the lowest process). 



We consider a system with periodic boundary condi- 
tions. This ensures translational invariance and the total 
momentum of the bound state will be a good quantum 
number. It is thus convenient to express the Hamilto- 
nian in the momentum basis \p, I) — e~~ lpl / 2 Y^j e~ lpj l), 
where p denotes the total momentum of the bound state 
(we use units where the lattice spacing is one) . Note that 
the ferromagnetic state necessarily has zero momentum 
\p = 0,0). In the momentum basis, the Hamiltonian is 
diagonal in p and acts as follows 

H\p,l) = (l-Si, ){[J t + h z l-2J ± cospS hl ]\p,l) 

- 2 J a cosp [|p, I + 2) + (1 - ^0(1 - 6 L2 )\p, 1-2)}} 
-J a (\p,2)S li0 + \p,0)S h2 )S pfi . (5) 

Because H a flips two spins, the state sectors with even 
and odd values of I are decoupled. Note that H± only 
affects the odd I sector. This follows from the fact that 
H± is only non-zero when it acts on the state with a 
single down-spin, the action on all other states produces 
more domain- walls, see Fig. [3] 

In order to solve the eigenvalue problem, we first 
parametrize the energy eigenvalues as 



E n (p) = J z + hzjJL n , 



(6) 



where p is the total momentum, n is a positive integer, 
which labels the internal excitation mode of the bound 
state, and fx n depends on p and coupling constants of 
the Hamiltonian. Because the excitation spectrum sepa- 
rates into distinct sectors, we will reserve the odd(even) 
n values for labeling the energy levels in the I odd(even) 
sector. The energy of the n'th mode is found by deter- 
mining the dimensionless quantity p, n from the following 
equations: 



(x) 



J_ 



z, n £ odd, 
0, n g even, 



(7) 
(8) 




FIG. 4: (Color online) The behavior of the ratio of Bessel 
functions J- V (x = 1)/Ji- V (x = 1) as a function of v is shown 
as the black solid curve. The blue dot-dashed line marks the 
value 7 = 2, while the red dashed curve shows the right hand 
side of Eq. ( |14[ | for y = 3. The red circle indicates the crossing 
point that gives the lowest energy solution vq — po/2 of the 
p — even sector. 



with J v {x) being the Bessel function of the first kind of 
order v. The upper(lower) equation is obtained by con- 
sidering the odd(even) I sector. Here, we have introduced 
the new variables 



2J a cosp Jj_ 
x = , z = —. 

h z J a 



(9) 



The even sector equation ([8]) is only valid for p ^ 0. The 
p = case will be considered separately below. 
Both equations are of the form 



J-u{x) 
Ji- V {x) 



7> 



(10) 



where 7 is a constant. Analytical solutions of this equa- 
tion has been found in some limits 16 . To get an intuitive 
picture of the solutions of Eq. (10 1, we have plotted the 



left hand side of the equation for a particular value of x 
in Fig. [4] as the black solid curve. For 7 — > 0, relevant 
for the even sector and the small z limit of the odd sec- 
tor, the solutions are gotten by the zero crossings of the 
curve. We see that they occur almost exactly at positive 
integer values of v, except for the lowest value which is 
somewhat below 1. In Fig. [5] we show these values of 
v as a function of x for 7 = 0. To a very good approx- 
imation, v is a positive integer as long as v > 1 + |x|. 
For the even I sector where /i = 2v, this implies that the 
solutions p, n , n 6 {2, 4, 6, . . .} are to a good approxima- 
tion even integers p, n — n when n > 2\x\ +2. For lower 
values of n, [i n is generally lower and depends on x. For 
small 7, relevant for the odd I sector where fi = 2v — 1 
when z is small, the solutions are very similar to the case 
7 = 0. This means that the solutions p, n , n £ {1, 3, 5, . . .} 
are odd integers /i„ — n for n > 2\x\ + 1. The qualita- 
tive effects on the solutions v of changing 7 for a fixed 
value of x can be inferred from Fig. |4j While the higher 
energy levels do not change substantially for this value 
of x, the lowest energy level decreases with increasing 
7. In Fig. |6j we have plotted the solutions of Eq. ( 10 ) 
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FIG. 5: (Color online) Solutions v of Eq. ( 10 1 as a function 



of x for 7 = 0. The red dashed line is the line v = 1 + \x\. 




FIG. 6: (Color online) Solutions v of Eq. |10| as a function of 
x for 7 = 2. The dot-dashed blue line has slope (7 + 7~ 1 )/2 
and the red dashed line is the line v = 1 + x for positive x. 



as functions of x for a fixed value of 7 = 2. For pos- 
itive x, the lowest energy state decreases with increas- 
ing x. This decrease becomes linear at large x and its 
slope can be found by writing the Bessel function ratio 
J- v (x)/Jx-v(x) = -[v/x+ ^lniJ-vix))]- 1 and using 
the asymptotic expression^ 



J_ v (-v6 



(11) 



valid for v — y 00 and < 8 < 1. We find that the lowest 
energy state solution decreases as v ~ — (7 + 7~ 1 ):e/2 for 
x — > 00 and 7 > 1. A line with this slope is overlaid 
on the plot in Fig. [6] for 7 = 2. The large x limit corre- 
sponds to the limit of vanishing magnetic field. Using the 
asymptotic result and the definitions of 7 and x, we get 
for the lowest excited state in the limit of vanishing h z , 



E = J z +h z (2v-l) ~ J z -2J ± (l + (J a /J ± ) 2 ) cosp corre- 
sponding to a spin wave excitation. This result can also 
be obtained by second order perturbation theory in the 
limit of vanishing magnetic field when J z 3> J± 3> J a . 

Extrapolating the spin- wave line to negative x, we see 
that it coincides smooth step-like behaviors of the energy 
levels. These steps become sharper and higher as 7 is 
increased and become a step discontinuity of unit one for 
7^00. This feature restricts the existence of the WZL 
for odd n to n > (7 + 7~ 1 )|a;| + 1 for x < and 7 > 1. 

From these considerations, it follows that a 
momentum-independent WZL 

E n = J Z + h z n, (12) 

with integer n, is restricted to high energies where 



J* 



^ 1 j_ 2J - 
n > 1 + - — x \ 

h z [J±/Ja + Ja/J±, 



J± < Ja, 
Jj_ > J a . 



(13) 



For p — 0, the sector with domain walls of even / will 
also couple to the ferromagnetic state. Thus, at p = 0, 
the equation for the energy levels in the even sector is 
different from the one at p ^ 0, Eq. For p = 0, we 
get 



J2{x ) 4(i/ + /x n )' 



(14) 



where y — J z /h z and xq = 2J a /h z . A similar equation 
including also the effects of an optical phonon at p = 
was obtained in Ref. 1131 The equation for p = has an 
additional negative solution /i < well separated from 
the other positive solutions. This solution arises from the 
singular behavior of the left hand side in the vicinity of 
|U = —y. Figure |4] shows a graphical solution of Eq. (14) 
with fj, — 2v for fixed values Xq = 1 and y — 3. Because 
of the rapid variations of the ratio of Bessel functions 
on the left hand side around positive even integer values 
of /i, we see that a relatively small finite value of the 
right hand side changes only slightly the even integer 
solutions found for p 7^ in the even sector. However, at 
low energies, there is a crucial difference. The right hand 
side has a singularity at fi = 2v = —y. The left hand 
side is positive for negative /i and increases as [i — > —00 
while the right hand side is negative for /1 > y before 
it changes sign as \x passes y. Thus somewhere below 
this singularity a negative solution fi will occur. This 
means that the ground state energy will be negative as 
J + h z fi — for fi = — y. For J a <C J z we find for the 
lowest energy solution approximately 



E = J z + h z no 



J z + 2h z 
6 



(15) 



x 1 



12J2 



■r 



{J z + 2h z y 



Jz + 2h z 



consistent with what is expected from second order per- 
turbation theory. 
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The energy eigenfunctions are given by \n,p) 
Y^hLo^n,i(p)\p,l} with coefficients 



lpn,l(p) 



!-(-!)< 



(16) 



valid for n odd, and 

•1 + (-!)' 



2 



[^a- M „)/2(a;)(l-^,o) 



£»(p = 0) 



^(2- A i„)/2(a ; ) ( 5p,o<5i,o 



(17) 



for n even. Note that only the odd(even) I coefficients are 
non-zero for odd(even) n. For small x, the Bessel func- 
tion is maximal when I = fi n w n. Qualitatively, this 
implies that for large magnetic fields, the n'th mode of 
the bound state is dominated by the state having n over- 
turned spins. This domination is total at p = 7r/2,37r/2 
where x = 0. 

These energy wave functions are orthogonal when the 
variable fi n obeys one of the Eqs. Q, pH), or (14 1 due to 
the Bessel function property^ 



E 

k=l 



J k +v(x)J k+fl (x) = 

X J\ +V {x) Jl +fi (x) ( J v [x) 



(18) 



v — /i 



Jy( x ) 

Ji +V {x) Ji+n(x) 



IV. DYNAMIC STRUCTURE FACTOR 

Having the eigenfunctions \n,p) and energies E n (p) we 
can calculate the dynamic structure factor, which at zero 
temperature is 

S aa \q,uo) = J2MS%\n,q)(n,q\S^' \0,0) 

n 

x5(u - (E n (q) - E (0)), (19) 

where |0,0) is the ground state and we have restricted 
the intermediate states to the states \n, q) in the zero and 
two domain-wall sectors that restricts the energy transfer 
uj < 2J Z . In the following we will consider S xx , S yy , and 
S zz separately. 

For S xx and S yy , it is convenient to consider the raising 
and lowering operators S ± — S x ± iS v . Expressed in 
terms of these the transverse dynamic structure factors, 
S xx and S yv are 



5" 



gyy 



[(s + - 



s-+) + (s- 
s-+)-(s- 



s ++ )] 
s ++ )] 



The ground state has zero momentum and can for J a <C 
J z be approximated by the ferromagnetic state \FM) 
where all spins point along the magnetic field. The 



calculation is simplified greatly by this approximation 
as then the structure factors S , S ++ , and are 
zero, which follows from S + \FM) — implying that 
gxx = S yy = 1$+-. The action of S~ on the ferro- 
magnetic state creates a state with one down-spin, thus 
it belongs to the odd sector, and will have rather high 
energy, of the order J z + h z . Using the eigenfunctions, 
we find 

oo 

S+-(q,oj)=J2^-E n (q))I n (q), 



where /„ is normalized relative intensity of the n'th 
mode, 



Uq) 



l^=l(g)l 2 



(20) 



Using the expression for the wave functions and the fol- 
lowing Bessel function identity 



l = l 



dv 



Ji -v-i{x) 

Jl -v{x) 



(21) 



with Iq an integer, the intensity can be expressed in the 
form 



d_ 



J l _ l /2+l(x) 



(22) 



where p n is the solution of Eq. ^ for the odd domain 
length I. 

For larger values of J a /Jz, it is no longer adequate 
to approximate the ground state with the ferromagnetic 
state. Taking into account the exact nature of the ground 
state gives additional contributions to S A and the cor- 
responding intensity becomes 



V i>0 / 

where we have used the following notation for the nor- 
malized wave functions: 



cm 



(24) 



We will omit the momentum label for the ground state 
Cf as it has zero momentum. The leading terms of the 
intensity give the contribution 

If (<?) « (C° ) 2 I n (q) + 4C° C° 2 C?(q)CS(q) cos q. (25) 

In this case we also get non-zero contributions to S 
and to S , S ++ which cause S xx to be different from 
S yy . In the same notation as above their contributions 
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FINITE TEMPERATURE 



4 ++ (<?) 



C?C "(0)5 3 , + 2e~W ]T C?Cf-i(9) cos(gZ/2) 
;>i 

C °C?(g) + 2^ C?CT +1 (?) coB(gi/2) j 



2>0 



x |c?C£(0)5 fl)0 + 2 £ cos(gJ/2) 
In-(l) = (# + (<z))*, 



where the * means complex conjugation. Approximating 
these with their leading terms, we get 

I-+(q) « (c?<tf(0)* g , + 2C§C?( g )cos?) 2 , (26) 
I++(g) « C?Q(0)(cgCr(0) + 2C°Cy(0))<y fl)0 



+2C C 2 °(C 1 i (g)f cosg. 



(27) 



Applying the operator S z to the ground state does not 
change the parity of I, thus all contributing intermedi- 
ate states have even n. It is convenient to split off the 
ground state contribution as it has zero momentum and 
frequency, it represents the squared magnetization, and 
write 



S"(q,u) 



6 qfi S(u) (28) 



where 



TO 



El 2 C?C?(q)sm(ql/2) 



sin 2 (g/2) 



(29) 



The leading contribution of the sum is 

4£ 2 |^= 2 (g)| 2 



4J„ 2 



cos 2 (g/2) 



(Jz + 2h z ) 2 



cos 2 ( g /2)C( 9 ), (30) 



where the introduced intensity J° v corresponds to contri- 
bution from the states with even n and can be written 



as 



TO 



(31) 



where fi n is the solution of the Eq. ([8]) for the even domain 
length I. 



At finite temperatures, the dynamic structure factor 
will in addition to transitions out of the ground state 
also get contributions that depend solely on the excited 
states. In particular, there will be contributions at low 
frequencies corresponding to the spacing between energy 
levels. It was these temperature induced contributions 
that were the focus of the neutron scattering experiments 
in Rcfs. 14151 We will consider such temperature-induced 
contributions at relatively low frequencies oj < J z . 

At finite temperature, entropic factors make it favor- 
able to break up a spin domain, thus inducing domain- 
walls. The magnetic field confines pairs of domain-walls 
leading to a picture of the finite temperature state as con- 
sisting of several bound states; short spin-down domains, 
interspaced by longer spin-up domains. To handle these 
thermal states, we will use the exact quantum mechanical 
treatment of an isolated bound state, and neglect the in- 
teraction between different bound states. We expect the 
quality of this noninteracting bound state approximation 
to be good on time scales shorter than the typical colli- 
sion time between bound states. This collision time can 
be estimated by the mean distance between the bound 
states, which is the typical size of a spin-up domain, di- 
vided by the velocity of a bound state. The typical length 
of a spin-up domain in units of the lattice spacing can 
be estimated from the emptiness formation probability 
for the Ising model in a magnetic fielcP^ and gives 



1 



1 — a-f ' 



where 



cosh(/3/i z /2) + y / sinh 2 (/3/i z /2) + e"^- 



(32) 



(33) 



where /3 is the inverse temperature. Using this and the 
maximum velocity of a spin down bound state v max , we 
expect the independent bound state approximation to be 
good for frequencies 

27TV max 



uj > 



(34) 



The bound state velocity v n (p) — dE Q^ pS> is largest 
for low-lying energy modes. For higher-energy modes, 
the dispersion becomes flatter and their velocity ap- 
proaches zero. For the n = 1 mode, the energy varies 
as Ei{p) w -2Jj_[l + (J a /J±) 2 } cosp for J± > J a at 
low momenta p, which implies a maximum velocity v\ = 
2Jj_[l + (J a /J±) 2 ]. For Jj_ < J a , the n = 1 mode be- 
haves almost as the n = 2 mode, which for 2J a /h z < 1 
has a maximum velocity t> 2 ~ 2J 2 //i z that increases for 
smaller fields and approaches t> 2 <~ 4J a . Thus the maxi- 
mum velocity of a bound state is v max = max(wi,w 2 ). 

The above validity criterion (34) takes into account 
the center of mass motion of the bound states. In ad- 
dition, the quantum mechanical uncertainty in the size 
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of a bound state can also ruin the noninteracting bound 
state approximation. From the wave functions, we esti- 
mate the size uncertainty to be ±2J a /h z , which implies 
that adjacent bound states have nonoverlapping bound- 
aries when 



2./,, 



4J a /h z < £ t . 



(35) 



Keeping in mind these restrictions, we can write down 
the dynamic structure factor at low frequencies U) < J z 
in the independent bound state approximation as 



few) 



xS[uj- {E m , ( P + q)- E m {p))] , (36) 



where 



S^' m {p,q) = (m,p\S« q \m , ,p + q)(m / ,p + q\S«'\m,p) 

and n m p is the occupation number of a bound state with 
internal energy index m and momentum p. Its functional 
form depends generally on the statistics of these excita- 
tions, but is expected to behave at low temperatures as 
n m>p « e"^ (Bm(p) " Bo(0)) K(^), where k is close to unity 
forV < J z . 

The action of the operator S z on a state with a spin 
down cluster of I spins and momentum p is 

S* \P, I) = y |p, 0^,0 - X T J ^ \P + Q,l), (37) 

which implies that the matrix elements S^ n is given by 
the expression 



s£„ = s q ,o ( y E c ?(p)Gr(p) - E ^ro)cr» 



+ (i-5, 



and the argument of the Bessel function is £ = smq| 
For q = 0, we get also a contribution from the ground 
state magnetization squared: 



This expression is rather difficult to deal with analytically 
for general values of the parameters. However, in the 
region of parameters where we expect the spectrum to 
be the WZL, we can evaluate it analytically. Focusing 
on this region where 2J a /h z <C 1 and Jj_ < J a , we get, 
see Appendix, 



S"(q,u) = 



N 



1 - e-"*- 



E G k (q)S(u-2h a k), 
(39) 



k=-N 



where the contribution from each mode for q =^ is 
J 2 (C) eP h - + l 



G (q) = 
Gk{q) = 



cosh(/3h z ) — cosq 2 

4(0 fi, ^>o, 

2sin 2 (g/2) le' 32 ' 1 ^, k < 0, 



G o (0) 
G fc (0) 



2Ja 



iV 


1 e' 3 ' 1 * + 1 


_ e -f>h z 


f 2 cosh(/% z ) - 1' 




fc = 1, 




fc = -1, 




|fc|>l. 



Thus there are contributions for uj — 2h z k, where k is 
an integer. The 2h z reflects the fact that S z does not 
change the parity of I. If we introduce the Bloch fre- 
quency lob — 2h z , this gives a result similar to the ex- 
pression found in Ref. [3] However, the temperature de- 
pendent factors are different. In particular, we get a prc- 



factor e 



which is a consequence of the occupation 



number of bound states. This is in contrast to the factor 
e -/3^ z /2 expected in the single domain- wall approxima- 
tion. A noteworthy feature of Eq. (39) is that increasing 



the magnetic field moves the spectral weight to lower 
Bloch frequencies. This follows from the fact that the 
maximum of the Bessel function squared Jk{C) f° r fixed 
£ occurs when k w £ — 1. 

For the transverse dynamic structure factor, we find 



Smn(P,<l) 



CZ(0)C?(q)d P;0 



+2EQ»Q+i(P + <7) cos 
;>o 



(40) 



ql-p 



This expression can also be evaluated analytically in the 
region where the spectrum is the WZL, see Appendix, 



S+- (q, u) = E%- h ^ k + !)) ( 41 ) 

k 

x2j 2 fn/ 1 ' fc>1/2 ' 

x ^ J fc^ j e /3M2fe+i) 5 k <\/2, 

where £ = 2 ^|sing| and k is an integer variable. For 
the transverse structure factor, the excitations occur at 
frequencies that are an odd multiple of h z , a consequence 
of the fact that S~ changes the parity of I, the number 
of overturned spins. 

VI. COBALT CHLORIDE 

C0CI2 • 2H2O is a quasi-one-dimensional anisotropic 
spin- 1/2 magnet, proposed in Ref. 3 as a candidate 
exhibiting BO in a magnetic field. C0CI2 • 2H2O 
has a dominant ferromagnetic coupling J z along the 
chains, which was determined from far-infrared absorp- 
tion spectroscopy to be J z — 36. 5K. Other intra- 
chain couplings J x and J y are smaller but non-zero. 



(38) 




FIG. 7: (Color online) Energy levels vs. momentum com- 
puted with the parameters in Eq. ( |42[ | for a magnetic field 
h z / Jz = 0.05. The dot-dashed blue curve corresponds to the 
asymptotic dot-dashed line drawn in Fig. ©but with 7 = 1.43, 
and the red dashed curve corresponds to the red dashed line 
in Fig. § 



The values of these couplings as well as other inter- 
chain couplings have been inferred both from far-infrared 
spectroscopy^ and fro m spi n wave analysis of neutron 
scattering experiment ;} ' ' I In this article, we use the 
following values to describe C0CI2 • 2H 2 0: 



J^=36.5K, J a = 3.8K, Jj_ = 5.43K. 



(42) 



An important consequence of interchain couplings in 
C0CI2 • 2H 2 is that they cause the spins to order anti- 
ferromagnetically below Tjv = 17.3 K. This implies that 
in the antiferromagnetic phase below Tjf the magnetic 
field h z used here should be interpreted as a sum of 
the external applied magnetic field and an internal field, 
which arises due to the magnetic moments of neighboring 
chains^. 

We have plotted the energy levels E n (p) for the above 
couplings in Fig. [7j The WZL is present at low energies 
in the momentum region around p = n/2 (3ir/2) and is 
bounded by the red dashed and blue dot-dashed curves, 
which correspond to the asymptotic lines drawn in Fig. [6] 
For energies E > J z + 2h z +4J a the spectrum is the WZL 
for p = and for E > J z + h z + 2J±[1 + (J a /J±) 2 ], it 
extends also to the region above ir so that the spectrum is 
the WZL for all momenta. For regions of energies where 
the spectrum is not WZL for all momenta, it is possible to 
see from Fig. [7] that the even levels are symmetric around 
7r/2 while the odd levels lack this symmetry property. 
This is a consequence of the skewness of levels seen in 
Fig. I] 

In order to see the effects of J± , we have in Fig. [8] also 
plotted the energy levels when J± — for comparisorP2 
We see that the main effect of J± is to lower the energy 
of the lowest odd level and to shift the low-energy odd 
levels in the region ir/2 < p < 3ir/2 so that they al- 
most coincide with the even levels. The even levels are 



FIG. 8: (Color online) Energy levels vs. momentum com- 
puted with J ± = 0, JajJz = 0.104, and h z /Jz = 0.05. The 
red dashed curve corresponds to the red dashed lines in Fig. [5] 
and marks the lower boundary of the WZL. 



unaffected by J±. 

The zero temperature transverse dynamic structure 
factor at h z /J z = 0.05 is shown in Fig. [9J Only 

transitions to odd n levels have nonzero intensity and it 
is seen that most of the spectral weight occurs for transi- 
tions to the spin- wave like state n = 1. The intensities of 
higher excited levels are weak in the momentum region 
around p = n/2 where we expect to see the WZL. Exactly 
at p — 7r/2, the size of the bound state is a good quan- 
tum number, thus higher excited bound states with 11 > 1 
have no amplitude to have the size 1 = 1, which is the 
dominant intermediate state generated by neutron scat- 
tering on a ferromagnetic state. Any intensity of n > 1 
levels at p = ir/2 reflects how the ground state devi- 
ates from being fully ferromagnetic. For the coupling 
constants relevant for C0CI2 • 2H 2 0, the probability for 
finding all spins up in the ground state is roughly 99%, so 
fluctuation corrections to the ground state are small and 
the integrated spectral intensity above the n = 1 mode 
is less than 1%. In Fig. 10 we show how the intensi- 



ties Eq. ( 23 1 of the different levels vary for two momenta 
q = and q = tt. For q = 0, the intensity drops ex- 
ponentially with frequency, but for q = tt, the intensity 
decreases only slightly before it increases up to the energy 
where the WZL sets in and then drops rapidly. This also 
reflects the fact that the main contribution comes from 
transitions to the I = 1 state. 

The behavior of the longitudinal dynamic structure 
factor S zz for parameters relevant for C0CI2 • 2H2O is 
shown in Fig. 11 Here, only excitations to even n levels 
are nonzero which implies that S zz is independent of Jj_. 
The total spectral weight of S zz (q 7^ 0) is however much 
smaller than for because it is proportional to the 

probability for finding two overturned spins in the ground 
state. This is reflected by the small factor J 2 /(J z +2h z ) 2 
in Eq. p0). 



For finite T, the validity of the noninteracting bound- 
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FIG. 9: Gray-scale plot of S + (q, uj) vs. q and u> at T = 
for J ± /J z = 0.149, J a /J z = 0.104, and h-JJ z = 0.05. 




0)/J z 



FIG. 10: Intensities 7 + "(g) vs. w at T = for Jx/J z = 
0.149, J a /J z = 0.104, and h z /J z = 0.05. The results for two 
momenta are shown, q = (open circles, dashed line) and 
q — it (solid circles, solid line). The lines are guides to the 
eye. 



0J/J z 



1.6 



1.4 





2?r- 



FIG. 11: Gray-scale plot of S zz (q, u) vs. q and lj at T = 
for J ± /J z = 0.149, J a /J z = 0.104, and /i z /J z = 0.05. The 
intensity of the plot has been increased by a factor 28 in order 
to make it visible on the same gray-scale as used in Fig. [9] 




FIG. 12: Region of validity of the noninteracting bound state 
approximation using the parameters in Eq. (421. The region 



where the inequality ( 35 1 holds is shaded in light gray (which 
also overlaps entirely the dark gray region). The region where 
the noninteracting approximation can be used for freq uen cies 
down to the Bloch frequency lob = 2/i z , inequality (341, is 
shown in dark gray. 



state approximation for C0CI2 • 2H 2 used here is con- 
strained mostly by J±. Its relatively large value causes 
the n = 1 bound state to have the largest velocity which 
according to the inequality (34) gives a lower bound on 
the frequency for which our approach is valid. If wc 
require that this lowest frequency equals the Bloch fre- 
quency lu b = 2h z , the noninteracting bound-state ap- 
proximation will be valid in the temperature/magnetic 
field region shaded dark gray in Fig.[l2]for the parameters 
relevant for C0CI2 • 2H2O. We do not expect our results 
to apply outside this region as a treatment of bound-state 
collisions is needed there. 

In order to see finite temperature signatures of BO, 
a high temperature is needed to thermally occupy the 
bound state levels. However, at high temperatures, the 
validity of our approach is restricted to large magnetic 

Increasing the magnetic field 



fields, as seen from Fig. 12 



has the disadvantage that the weights of the finite fre- 
quency Bloch peaks become small Eq. (39), thus weak- 



ening the signatures of BO. Therefore a judicial choice of 
temperature and magnetic field must be made to make 
observations possible. 

The optimal magnetic field for the first resonance at 
ujb is h z ~ J a . We will use a larger magnetic field, 
h z = 0.2 J z , as that allows our approach to be used up 
to a temperature T « J z /2. We find that for S zz the 
maximum intensity of the finite frequency WZL transi- 
tions occur at q = ir/2. In Fig. 13 we have plotted our 



analytical result Eq. ( 38 1 numerically. We compare this 



with the expression obtained in the WZL limit Eq. ( 39 ) 



using the same parameters. We see a clear peak at the 
Bloch frequency ug — 2h z also when the conditions for 
the WZL are suboptimal as is the case with the parame- 
ters in Eq. ( 42 ) . The WZL calculation (red dashed line) 



overestimates the weight of the peaks, but do reasonably 
capture their relative intensities. For higher tempera- 
tures, the WZL expression matches Eq. ( 38 1 better as 
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S zz (n/2) 




FIG. 13: (Color online) S zz (q = ix/2,ui) vs. lo/J z for T = 
J z /2, h z /J z = 0.2, and parameters from Eq. (421. The red 
dashed curve is the corresponding WZL result Eq. ( |39| using 
the same parameters. In order to generate the plot delta- 
functions were approximated by a Gaussian distributions with 
variance 10" 4 . The vertical axis values are in multiples of 
k{P = 2/Jz), a number of order unity. 



S + -(n/2) 



12 






10 


[ i 




8 
6 


i 

\ i 




4 




ii 

i 1 1 


'« 2 




i 1 1 
i i i 



-0.4 -0.2 



0.2 0.4 0.6 0. 



0J/J z 



7t/2,cj) vs. uj/J z using 
The red dashed curve is 



FIG. 14: (Color online) S + ~(q = 
the same parameters as in Fig. |13| 
the corresponding WZL result Eq. ( 41 1 using the same pa- 
rameters. Delta-functions were approximated by a Gaussian 
distributions with variance 10" 4 . The vertical axis values are 
in multiples of k(/3 = 2/J z ), a number of order unity. 



then more emphasis is put on the higher-energy part of 
the spectrum which is more WZL-like for all momenta. 
We wish to emphasize that the thermally induced transi- 
tions here come with an overall factor e~P J * which makes 
them difficult to observe at low temperatures. 



In Fig. 14 we have plotted S + using the same param- 



eters as in Fig. 13 Peaks at frequencies corresponding to 
odd multiples of the magnetic field arc clearly seen among 
other peaks caused by the dispersion of the lowest-energy 
modes. 



VII. CONCLUSION 

In this work, we have investigated the possibility 
of observing spectral signatures of magnetic BO in a 
one-dimensional anisotropic ferromagnetic spin system 
placed in a magnetic field. This system was considered 
previously^ but within an approximation where only a 
single domain-wall was included. We argue that the sin- 
gle domain-wall approximation is insufficient at a finite 
magnetic field. Instead, we consider a bound state of a 
domain-wall and an anti-domain wall; a spin cluster of 
adjacent spins antialigned with the magnetic fielcP, and 
treat the thermal state as a noninteracting gas of such 
excitations. This allows us to also include the effects of 
the additional coupling Jj_, which probably is present in 
most anisotropic materials that have a nonzero value of 
J a . For instance in C0CI2 • 2H2O for which neutron scat- 
tering searches for BO have been made, J± is bigger than 

Ja- 

We have treated the quantum mechanics of the bound 
state and obtained its energy levels and wave functions. 
The spectrum in a magnetic field will be split by the 
magnetic field essentially into modes corresponding to 
the size (number of overturned spins) of the bound state. 
In the momentum region around p = tt/2 (p = 3n/2), 
the energy levels are equidistantly spaced down to the 
lowest energies. This corresponds to the WZL. For other 
momenta the exchange couplings, J a and J±, cause dis- 
persion of the low-lying energy levels. This effect dimin- 
ishes for higher energies and the spectrum becomes the 
WZL above a threshold energy for all momenta. 

We have also calculated the neutron scattering dy- 
namic structure factor at zero and low temperatures. At 
zero temperature, the expected response occurs at high 
frequencies and there should be considerable chances of 
seeing the magnetic field splitting of the spectrum. How- 
ever, it will be difficult to see the WZL because of the low 
spectral weight in this region. This is because neutron 
scattering flips a single spin and couples most strongly 
to the bound states that have a significant amplitude of 
having a single down-spin cluster, which has a disper- 
sion that is heavily influenced by J±. The longitudinal 
channel is not influenced by J±, but is much weaker at 
non-zero momentum. 

Detecting BO at finite temperatures with neutron scat- 
tering seems more promising. However, the thermal 
occupation number of bound states implies that these 
signatures will be suppressed at low temperatures as 
e —pJ* _ This might be the reason for the nonobservation 
of BO in the neutron scattering experiments^'. Never- 
theless, it might still be possible that neutron scattering 
on C0CI2 ■ 2H2O can be used to observe signatures of 
BO provided a careful selection of temperature and mag- 
netic field is being made. In that respect, care must be 
taken so as to secure a big enough thermal population of 
bound states, a large intensity of the finite frequency res- 
onance^), and a regime where collisions of bound states 
do not destroy the BO. 



11 



Appendix A: WZL structure factors 

Let us calculate the dynamic structure factors at finite 
temperature in the limit where the spectrum is the WZL, 
that is h z ^> 23 a and J± <C J a - In this limit E n (p) = 
Jz + h z n, and the wave functions are Bessel functions of 
integer order 



cT(p) = 



! + (-!) 



l-n 



- J i- n (xq cos p) , 



(Al) 



where xq = 2J a /h z <C 1. In the limit xq these 
wave functions are normalized to unity. This also holds 
approximately at small finite xq which follows from the 
property of Bessel functions: Y^k=- oo 4( x ) = 1 anc ^ 
the fact that contributions from higher Jk{xo) decrease 
rapidly with fc for small xq. Thus only small errors are 
introduced by extending the sum to — oo at small finite 
x . 



For S zz (q ^ 0,w), the matrix elements (38) become 



Wo sm 2 (q/2) 



J i- n (xq cosp) J i-m (xq cos(p + q)) sin 



z>o 



(A2) 



ql 



where the sum over I goes over even(odd) integers when 
n and to both are even(odd), otherwise every term in the 
sum is zero. It is convenient to introduce a new integer- 
valued variable t — (I — n)/2 and rewrite the sum in the 
form 



where integration over momentum p was performed for 
integer values of n and to variables 



dp 

— sn 
2tt 



q 2 



1/2, 



n 7^ m, 



sin (qf ), n — to. 



No transitions between the even and odd sectors are al- 
lowed, thus it is convenient to introduce a new integer- 
valued variable fc = (m — n)/2 which describes the energy 
difference between the states involved in the transition. 
We can reorder the double sum as 



e = m^m-m: 

m,n rn>n m<n n rn>ri n 

+EE =£E+EE> 



(A6) 



n k>0 m k<0 



that allows us to rewrite the dynamic structure factor 
in the case of non-zero energy transitions, k ^ 0, in the 
following form 



S* 



1% 2sin 2 (g/2) e? h - 



(A7) 



mo r 



k > o, 

fc < 0, 



where we used the expression for the sum of the first N 
terms of a geometric series 



A? 

E> 



_ -I3h, 



1 _ e -M 



1 



o-Ph z 



1 



(A8) 



Jf(xncosp) J t | n-m (xn cos(jp+q)) sin 



t>-3 



Jn-mCC) Sin 



7T \ n — TO 

P ~ 2J 2 



q 2 



2 

, (A3) 



where we introduced a new variable C — x a \ sinq|. The 
sum over the product of Bessel functions was performed 
by extending the sum to negative — oo which only induces 
small errors when xq <C 1, and then using Graf's addition 
theorem^ 

oo 

J2 Jk+u(u) J k (v) ^(fcfl = JAw) ^>x) (A4) 

fc — — oo 



with the relations w = a/ it 2 + v 2 — 2uv cos 0, w cos x = 
u — v cos 4>, and iu sin x = v sm </>■ 

The dynamic structure factor S zz for q ^ in the inde- 
pendent bound state approximation can then be written 
in the form 



S zz {qM 



«C8) 



<z#o shr(g/2) 



5^ e -«J«+A-'») ( y( w _ ( m -n)) 



xJ^(C) 1/2 2 ' , (A5) 
2 sm (gfj, n — m, 



For the zero mode, fc = 0, in order to find the sum in 
Eq. (|A5| we can use the following identity 



oo -i*2 

SixT (an)e bn = r 



1 + e 



-b 



2 cosh — cos 2a 

n— l 

(A9) 

which can be proved using Euler's formula and sum of 
terms of geometric series. This gives the contribution to 
the zero mode 



1 _|_ P -Ph z 

lt° l-e-l 311 * 2 cosh(J3h z ) - cosg' 



8(u>) 



J 2 o(0 



(A10) 

Finally, combining together Eqs. (A7) and ( |A10 | we ob- 
tain the dynamic structure factor: 



S"{q,u>) 



K{p)e 



-/3(Jz+h z ) 



N 



1 



-0h z 



J2 G k (q)8(uo-2h z k), 
(All) 



k=-N 



where the contribution from each mode is 
J 2 (C) e? h *+l 



Go = 
G k = 



cosh(/3h z ) — cos q 

4(0 fi, 



fc > o, 



2sin 2 (g/2) |e^ fe , fc < 0, 



(A12) 
(A13) 
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and the argument of the Bessel function is £ = sinoj. 

The leading contribution to the dynamic structure 
factor comes from transition between states with 

nonzero momentum p. Then the matrix element in 
Eq. (40) becomes 



Since the only allowed transitions are between different 
parity sectors it is convenient to introduce the integer- 
valued quantity k = (to — n — l)/2. The integral over 
momentum is 



(A14) 



x J i+i- m (xq cos(p + q)) 



ql-p 



where the sum over I is over even(odd) integers when n is 
even(odd) and to is odd(even). Introducing the new inte- 
ger variable t = (I — n)/2 and summing over the product 
of Bessel functions using Graf's addition theorem gives 

^2 J f (xn cosp) J t , n-m+i (ffn cosfa + q)) 



t>-^ 



= J,i-r. 



- (C) cos 



-f] 

7r \ n — m + 1 p 



2 q 2 



dp 2 
— cos 
o 2tt 



7r \ p n 
2- p ) k+ 2 +q 2 



After reordering of the double sum we obtain finally 



S+-(q,oj) 



K(/9)e- 



-PJ. 



aPh z 



— J2s(uj~h z {2k+l)) (A16) 



x2j 2 fnl 1 ' fc>1/2 ' 



(^-15) where ( = \ s'mq\ and k is an integer. 



1 E. E. Mendez, F. Agullo-Rueda, and J. M. Hong, Phys. 
Rev. Lett. 60, 2426 (1988). 

2 M. Ben Dahan, E. Peik, J. Reichel, Y. Castin, and C. 
Salomon, Phys. Rev. Lett. 76, 4508 (1996). 

3 Jordan Kyriakidis and Daniel Loss, Phys. Rev. B 58, 5568 
(1998). 

4 W. Montfrooij, G. E. Granroth, D. G. Mandrus, and S. E. 
Nagler, Phys. Rev. B 64, 134426 (2001). 

5 N.B. Christensen, K. Lefmann, I. Johannsen, and O. 
J0rgensen, Physica B 276, 784 (2000). 

6 Z. Cai, L. Wang, X. C. Xie, U. Schollwock, X. R. Wang, M. 
Di Ventra, and Y. Wang, Phys. Rev. B 83, 155119 (2011). 

7 R. Coldea, D. A. Tennant, E. M. Wheeler, E. Wawrzynska, 
D. Prabhakaran, M. Telling, K. Habicht, P. Smeibidl, and 
K. Kiefer, Science 327, 177 (2010). 

8 Jonas A. Kjall, Frank Pollmann, and Joel E. Moore, Phys. 
Rev. B 83, 020407(R) (2011). 

9 J. B. Torrance Jr. and M. Tinkham, Phys. Rev. 187, 587 
(1969), Phys. Rev. 187, 595 (1969). 

10 J. Villain, Physica B+C 79, 1 (1975). 

11 H. Yoshizawa, K. Hirakawa, S. K. Satija and G. Shirane, 



Phys. Rev. B 23, 2298 (1981). 

12 S. E. Nagler, W. J. L. Buyers, R. L. Armstrong, and B. 
Briat, Phys. Rev. Lett. 49, 590 (1982); J. Magn. Magn. 
Mater. 31, 1213 (1983). 

13 H. C. Fogedby, Phys. Rev. B 10, 4000 (1974). 

14 S. E. Nagler,W. J. L. Buyers, R. L. Armstrong and B. 
Briat, Phys. Rev. B 28, 3873 (1983). 

15 N. Ishimura and H. Shiba, Prog. Theor. Phys. 63, 743 
(1980). 

16 S. B. Rutkevich, J. Stat. Mech. (2010) P07015. 

17 M. Abramowitz and LA. Stegun, Handbook of Mathemat- 
ical Functions, (Dover, New York, 1972). 

18 A. P. Prudnikov, Yu. A. Brychkov, and O. I. Marichev, 
Integrals and Series, 2: Special Functions (Gordon and 
Breach, New York, 1990). 

19 M. Suzuki, Phys. Lett. A 301, 398 (2002). 

20 J. K. Kjems, J. Als-Nielsen and Hans Fogedby, Phys. 
Rev. B 12, 5190 (1975). 

21 H. Shiba, Prog. Theor. Phys. 64, 466 (1980). 



